Late Reheating, Hadronic Jets and Baryogenesis 



Takehiko Asaka, 1 Dmitri Grigoriev, 2,3 Vadim Kuzmin, 2 and Mikhail Shaposhnikov 1 

institute of Theoretical Physics, Swiss Federal Institute of Technology, CH-1015 Lausanne, Switzerland 
Institute for Nuclear Research of RAS, Moscow 117312, Russia 
3 Mathematical Physics, Natl. Univ. of Ireland Maynooth, Maynooth, Co. Kildare, Ireland 

(Dated: October 8, 2003) 

If infiaton couples very weakly to ordinary matter the reheating temperature of the universe can 
be lower than the electroweak scale. In this letter we show that the late reheating occurs in a highly 
non-uniform way, within narrow areas along the jets produced by ordinary particles originated from 
infiaton decays. Depending on infiaton mass and decay constant, the initial temperature inside the 
lumps of the overheated plasma may be large enough to trigger the unsuppressed sphaleron processes 
with baryon number non-conservation, allowing for efficient local electroweak baryogenesis. 
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Introduction. — The inflationary paradigm (for a re- 
view see, e.g. books p]) happened to be very successful 
for understanding of the basic properties of the universe. 
It is assumed that the energy density of the early uni- 
verse was dominated by a potential energy of a scalar 
field - infiaton. The accelerating expansion of the uni- 
verse then leads to a solution of the horizon and flatness 
problems, whereas quantum fluctuations of the infiaton 
field result in density perturbations necessary for struc- 
ture formation and seen as the temperature fluctuations 
of the cosmic microwave background radiation. 

The exponential expansion of the universe during in- 
flation must be eventually replaced by a radiation dom- 
inated epoch, that started at least somewhat before nu- 
cleosynthesis and lasted till about recombination. The 
process of transfer of the energy density of the infiaton to 
ordinary matter is usually called reheating, and several 
scenarios for how it may proceed have been proposed. 
Quantitatively, they can be distinguished by the value 
of the reheating temperature Tr, below which the uni- 
verse expansion is dominated by radiation. If the cou- 
pling of infiaton to the ordinary matter is sufficiently 
strong, the energy transfer occurs very rapidly due to the 
phenomenon of broad parametric resonance right after 
inflationary stage Q; this leads to high reheating tem- 
peratures Tr ~ 10 10 GeV or so. If, on the contrary, 
the coupling is very weak, the exponential expansion of 
the universe is first replaced by a matter dominated pe- 
riod during which infiaton oscillates without dissipation. 
Then, the perturbative decays of infiaton heat the uni- 
verse up to some temperature Tr. In the latter case the 
reheating temperature can be very low, with the only re- 
liable bound Tr > 1 MeV coming from the successful pre- 
dictions of the big bang nucleosynthesis. The small values 
of the reheating temperatures can naturally occur in cer- 
tain supergravity models (see, e.g. where the decay 
rate of the infiaton is suppressed by the Planck scale. 

A successful cosmological model should also explain 
the absence of antimatter in the universe and the baryon 
to entropy ratio ur/s ~ 9 • 10 -11 0]. The baryogene- 



sis must occur after inflation since otherwise all created 
baryonic excess will be exponentially diluted. 

Quite a number of different baryogenesis mechanisms 
exist in theories with high reheating temperatures. The 
explosive particle production during infiaton decay in 
the wide resonance case p| may lead to production 
of GUT leptoquarks, the subsequent CP-violating and 
baryon number non-conserving decay of which gives rise 
to baryon asymmetry of the universe. The gravitational 
production of superheavy particles has the similar effect 
Q. The thermal or non-thermal || production of 
heavy Majorana neutrinos may prepare suitable initial 
conditions for generating the lepton asymmetry |9j , which 
is transformed then to baryon asymmetry due to anoma- 
lous electroweak number non-conservation |icj . If the 
reheating temperature is greater than the electroweak 
scale an electroweak baryogenesis (for reviews see jjjj) 
can take place. 

The theories with high Tr may, however, be in conflict 
with observations because of overproduction of danger- 
ous relics like gravitinos From this point of view 
the theories with small reheating temperature are more 
advantageous, as the production of unwanted particles is 
automatically suppressed. At the same time, the prob- 
lem of baryon asymmetry of the universe in much more 
difficult if Tr is relatively small, simply because in this 
case the low energy baryon number nonconservation is 
required. Thus, none of the mechanisms related to the 
production of heavy leptoquarks or Majorana neutrinos 
is operative, and one has to rely on some variant of elec- 
troweak 0, Q or Affleck-Dine baryogenesis. 

The electroweak baryogenesis, occurring in expanding 
and almost equilibrium plasma is highly constrained, as it 
requires the freezing of the sphaleron processes after the 
first order phase transition. This condition can be con- 
verted in the upper bound on the Higgs mass in the mini- 
mal standard model 0] . This bound cannot be satisfied 
with the experimental value of the top quark mass |l7j . 
so that new physics is required. In the minimal super- 
symmetric standard model there exist a (small) region of 
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the parameter-space leading to a sufficiently strong first 
order phase transition (see |18| and references therein). 

In [14J it was pointed out that the maximum temper- 
ature during reheating can in fact be much larger than 
the temperature Tr and that the rate of the universe ex- 
pansion in the region of the electroweak phase transition 
can be faster than it is usually assumed. This allowed 
the authors to relax somewhat the Higgs mass bound for 
the electroweak baryogenesis. 

In this letter we will show that the late inflaton de- 
cays heat up the plasma in a very non-uniform way and 
that the local temperature along the trajectories of de- 
cay products is substantially higher than the average one. 
For a range of parameters the mean temperature of the 
plasma is small enough to shut off the sphaleron transi- 
tions whereas the temperature of the overheated regions 
is large enough to switch them on. Electroweak baryo- 
genesis is then possible in these regions. The overheated 
regions cool down due to diffusion and expansion. This 
process has a highly non-equilibrium character, so no first 
order phase transition is required to satisfy the wash out 
condition and thus no bound on the Higgs mass is im- 
plied at all. Remarkably, the resulting baryon asymmetry 
does not depend much on many details of the process and 
may be consistent with the observed one for large enough 
CP-violation. 

Local overheating. — Let is the inflaton mass and 
= ftpM^ is its width. Assuming the instantaneous 
decay of the inflaton at time — l/r^ the reheat- 
ing temperature is given by Tr — -J J^M^Mq, where 
i 

Mo = Mpi/1.66gi , Mpi is the Planck mass, g* ~ 100 
is the effective number of massless degrees of freedom. 
The typical numbers we will be interested in are Tr < 
Tw ~ 100 GeV (here Ty/ is the freezing temperature of 
the sphaleron processes), and > 10 10 GeV, what re- 
quires rather weak decay constant ~ 10" 24 , or smaller. 
We shall assume, for simplicity, that inflaton decays into 
quark- antiquark pair, though other decay channels lead 
essentially to the same result. 

The number density of inflatons decreases with time as 
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?ioexp(— Tj,t) with scale factor a. The first de- 
cays occur essentially in the vacuum, whereas at t ~ t& 
inflatons are surrounded by the plasma with the temper- 
ature T ~ Tr. The decay products of inflaton, ultra- 
relativistic quarks with the energy M^/2 3> Tr, are in- 
jected into the plasma and heat it locally. Our aim now 
is to understand the typical size and geometry of the 
overheated regions, as well as their temperature. 

The dynamics of energy losses of high energy quarks 
and gluons in quark-gluon plasma is a complicated prob- 
lem which must incorporate the Landau-Pomeranchuk- 
Migdal effect 0] and non-abelian character of interac- 
tion of quarks and gluons. The main features of it have 
been understood only recently (see |20| and references 
therein), the crucial one being that the energy loss per 



unit length is increasing in infinite plasma as a square 
root of parton energy, see eq. below. 

The most significant part of energy losses is related to 
the soft gluon emission in the multiple scattering of the 
hard parton on the particles of the plasma. The energy 
spectrum I of the emitted gluons per unit length has the 
following approximate form for lu bh <C < Eg : 
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where u and E denote the gluon and parton energies, 
a s is QCD coupling constant, X g is the mean free path 
of the gluon, and o; BH is the Bcthe-Heitler frequency, 
w BH = X g n 2 . Here [i is the typical screening mass which 
is assumed to be of the order of the Debye mass in the 
plasma. From one gets the stopping distance of the 
initial parton L tot : 
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where 



w B h log (E o/o> BH ) for E > uj bb . The aver- 



age energy of the emitted gluons is (u>) = \/lu B¥1 Eq/2 and 
their number is given by N g = 2^/Eq] ^BH ■ 

If the emitted gluons have energies to 3> w BH , they lose 
their energies as the parent parton does. This cascade 
process will terminate when the energy of the emitted 
gluons becomes comparable to or smaller than uj bii . We 
find the energy for the n-th gluons as 
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so number Nrh of the cascade steps from (uj) = Eo to 
Mat bh < w B „ is typically 3-4 for E - 10 10 -10 n GeV. 

To find the geometry and volume of the region where 
the emitted gluons deposit their energy we note that the 
radiation with the frequency u> is mainly concentrated in 
the cone with a small angle 9, which can be computed 
with the help of the results of [21| and is given by 
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where is length of the energy loss for the gluon with 
energy u, 
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lo bh = lu bb Iyl(u)/lu bh ). The transverse distance r w = 
LujOu, traveled by a gluon with the energy uj is then 
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(the log factor should be neglected at u> ~ oj bh ). It is seen 
that, although the cascade process produces gluons hav- 
ing various energies, the largest r u is due to the gluons 
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with the lowest energy uj ~ a; BH - Thus, the overheated re- 
gion has an approximate cylindrical form with the length 
given by eq. J3J) and with the radius r* ~ r UBH . Its vol- 
ume is V = nrlL tot . 

Now we are at the point to estimate an effective tem- 
perature inside the cylinder in which the emitted gluons 
are living. By assuming the energy conservation and the 
rapid thermalization, we obtain that the effective tem- 
perature T», found from ^gfpT* = -f^, is given by 



T* ~ 5.3 x 10" 
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For numerical estimates we take the non-perturbative 
value of the Debye screening mass found in J^. As for 
the gluon mean free path, we take the value found in j2^| 
and denoted as A°' or simply take a gluon damping rate 
j g from [24|]. They are different by roughly one order 
of magnitude, what allows one to get an estimate of the 
uncertainties in Q. 

For example, when Tr — Mz and Eq — 10 11 GeV, 
we find (cf. Fig. ^| that the effective temperature is 
T* ~ 212 GeV (A s = l/ lg ) or T* ~ 151 GeV (A fl = A; 1 )- 
Fig- n clearly demonstrates that the local temperature of 
the overheated regions can substantially exceed the freez- 
ing temperature of the sphaleron processes, provided the 
background cosmic temperature is relatively high, say 
Tr ^ 10 GeV and the energy of the parent parton is ex- 
tremely high, say Eq ^ 10 10 GeV. 

Baryogenesis. — Let us assume now that the parame- 
ters of the inflaton are such that the temperature T* is 
high enough: T* > T aph - 100 GeV > T w . Here T aph is 
the temperature above which the rate of the sphaleron 
transitions is unsuppressed and is given by (per unite 
time and unite volume) r aph ~ na^yT^, where k ~ 10 
[25|. At the same time, the reheating temperature Tr 
can be small enough, so that baryon number is conserved 
away from the overheated regions (i.e. Tr < Tw). This 
highly non-equilibrium situation is possible at any choice 
of parameters of the underlying electroweak theory, and, 
therefore, the wash out bound of 0] is not applicable. 

The baryon asymmetry of the universe can be esti- 
mated as the number of sphaleron transitions which takes 
place inside the overheated regions and go asymmetri- 
cally due to CP-violation: 
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where n par t on is the number density of the highly ener- 
getic partons coming from the inflaton decay, At denotes 
the lasting time of the rapid shaleron processes and we 
have introduced <5 CP to represent the effective magnitude 
of CP-violation. 

With the two-particle decays of the inflaton the num- 
ber of partons is simply rt parton = 2rt0, and ^f- ~ 
from the condition defining the reheating temperature. 
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FIG. 1: Contour plots of effective temperature T* in the Tr- 
Eq plane for A 9 = l/"f g (up) and X g = Xg (down). 



Putting all factors together, we get 
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Quite amazingly, besides the expected CP-violating fac- 
tor, the result depends just on the reheating temperature 
and on the sphaleron transition time. In particular, the 
temperature of the overheated regions has canceled out 
from eq. What is important is that the overheated 
regions must be in the symmetric phase of the elec- 
troweak theory, where baryon number non-conservation 
is not suppressed. 

The overheated regions cool down by the growth of 
the volume due to diffusion, and the rapid spharelon pro- 
cesses terminate eventually when the temperature inside 
becomes T ph . From the energy conservation the radius of 
the overheated region at T Bph is given by I — r*(T*/T aph ) 2 
and is reached after the diffusion time At ~ £ 2 /AD (we 
find from j2|| that the diffusion coefficient D ~ l/7 9 ~ 
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O.lAjj 1 ), which gives At ~ £^ y^h) ' ^ or a re §i° n °f 
parameters At is long enough to thermalize VU-bosons 
which are essential to the sphaleron processes. 
Finally, one gets 

^>(1-10)-10- 7 «S CP , (10) 
s ~ 

depending on the estimate of the gluon mean-free path 
discussed above. So, the extension of the standard model 
with suitable CP-violation may work. 

Conclusions. — We have shown that a successful elec- 
troweak baryogenesis can take place in inflation models 
with low reheating temperature. Though our estimates 
are rather rude (even parton losses in plasma have rather 
large uncertainties due to the lack of exact knowledge of 
the kinetic coefficients), it is clear that the increase of the 
inflaton mass makes the temperature of the overheated 
regions higher, what triggers the mechanism at some crit- 
ical inflaton mass. 

Generally speaking, the plasma overheating by the pro- 
cesses discussed above will take place along the trajecto- 
ries of decay products of any sufficiently heavy particles 
provided their decay products interact with the back- 
ground plasma. 
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